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ABSTRACT 

We attempt to measure the proper motions of two magnetars — the soft gamma-ray repeater 
SGR 1900+14 and the anomalous X-ray pulsar IE 2259+586 — using two epochs of Chandra observa- 
tions separated by ~ 5 yr. We perform extensive tests using these data, archival data, and simulations 
to verify the accuracy of our measurements and understand their limitations. We find 90% upper 
limits on the proper motions of 54masyr~ 1 (SGR 1900+14) and 65masyr _1 (IE 2259+586), with 
the limits largely determined by the accuracy with which we could register the two epochs of data and 
by the inherent uncertainties on two-point proper motions. We translate the proper motions limits 
into limits on the transverse velocity using distances, and find v± < 1300 kms -1 (SGR 1900+14, for 
a distance of 5kpc) and v± < 930 km s -1 (IE 2259+586, for a distance of 3kpc) at 90% confidence; 
the range of possible distances for these objects makes a wide range of velocities possible, but it seems 
that the magnetars do not have uniformly high space velocities of > 3000 km s -1 . Unfortunately, our 
proper motions also cannot significantly constrain the previously proposed origins of these objects in 
nearby supernova remnants or star clusters, limited as much by our ignorance of ages as by our proper 
motions. 

Subject headings: astrometry — stars: pulsars: general — stars: individual: alphanumeric: 
(SGR 1900+14, 1E2259+58.6) — stars: neutron — X-rays: stars 



1. INTRODUCTION 

The sources know n as magnetars (see 
IWoods fc~T hompson 2006 for a review) , comprised of the 
observational classes soft gamma-ray repeaters (SGRs) 
and anomalous X-ray pulsars (AXPs) , are thought to be 
young neutron stars w i th magnetic field s abov e 10 14 G 
([Duncan k Thompson! Il992t IPaczvhskil 1 19921 ) ■ Why 
some neutron stars are magnetars (having magnetic 
field decay as their primary energy source and usually 
emitting no detectable radio pulsations) while others 
are ordinary rotation-powered rad io pulsars (with com- 
parat ively little X-ray emission; IKaspi feM cLaughlin 
120051 and references therein) is a fundamentally unre- 
solved issue, especially since some radio pulsars have 
been discovered with fiel d strengths comparable to 
those of magnetars ( e .g.. IKaspi k McLaughlin! 120051 
also see iGavriil et all 120081 ). The proposed expla- 
nation s run the g a mut fr om evolutionary sequ ences 
(e.g., iLin k Zhang 120041: iKasumov et al.l I2006 D and 
quiescent states (e.g., I Woods k Thompson! [20061 : 
iGotthelf et al.l I2004D to differences in progenitor mass 
(|Gaensler et all I2005D or magnetic f i eld orientation 
and geo metry ( | Zhang k Harding] l2000t iKulkarni et al.l 
2003). [D uncan k Thompson! ([19921 ) have proposed 
that an initial short birth period is responsible for 
the generation of the high magnetar fields through 
an efficient large scale dynamo, and have also pointed 
out that high magnetic fields can produce asymmetric 
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neutrino emission at birth, resul ting in extreme space 
veloci ties ( > 10 3 kms _1 ; also see iThompson k Duncan! 
119931 119951: lLalll2001[ ). Such large velocities would also 
help explain the apparent offsets between some SGRs 
and their putative n atal supernova remnants (e.g., 
Rot hschild et al.l Il994f) . These results for magnetars, 
taken with the exte nsive investigation of rotation- 
powered pulsars (e.g.. iNg fc Roman! 120071 ) . imply that 
the initial spin periods, surface magnetic fields, and 
birth kick velocities of NS all originate in supernova core 
collapse processes, and the physics of these phenomena 
are tightly interwoven. Magnetars, by virtue of their ex- 
treme magnetic fields, may provide a direct probe of this 
interdependence. For example, if birth kicks are driven 
by asymmetric neutrino emission mediated by high 
B- fields (a mechanism that may not b e able to produce 
the highest kick velocities; lLa.il l200lh . then magnetars 
should have a population velocity much higher than that 
of the radio pulsar population. However, a more critical 
investigation of some o f the claimed SGR/ supernova 
remnant associations led I Gaensler et al.1 ([200 1[ ) to doubt 
many associations and to conclude that magnetars as 
a class had velocities < 500kms _1 , comparable to the 
radio pulsar population. 

Such speculations were largely without concrete ex- 
amination, as the velocities of magnetars are difficult 
to measure. As a class, they lack radio counterparts 
and so cannot be used for traditional Very Long Base- 
line Interferometry, and their spin-down properties are 
too noisy for measurement of "timing" proper motions. 
This has changed recently, though. The first measure- 
ment of the direct proper motion (and hence space ve- 
locity, with an assume d dista nce) of a magnetar was the 
work of iHelfand et al.l (|2007D . who found a proper mo- 
tion of 13.5 ± l.Omasyr -1 (« 212 kms -1 at a distance 
of 3.5 kpc) for XTE J1810-197. This was made possible 
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by the detection of a bright, transient radio counterpart 5 
(jHalpern et al.|[200l . 

In this paper, we attempt to further address this issue 
by trying to measure the proper motions of two magne- 
tars, SGR 1900+14 and IE 2259+586. In the absence of 
steady or widespread radio counterparts, we try to mea- 
sure proper motions in the X-rays: while IE 2259 +586 
has an infrared counterpart (|Hulleman et al.ll2001l ) and 
would also be amenable to a measurement using data 
at those wavelengths, we try to use a technique appli- 
cable to objects without infrared counterparts, although 
we anticip ate infrared data being used to check our mea- 
surement ([Cameron. Britton. fc Kulkarni 20QM). 

Not only can our data provide valuable constraints 
on the velocity distribution of magnetars, but they can 
also potentially address associations between the mag- 
netars and other objects (supernova remnants and mas- 
sive star clusters). Such associations are valuable for the 
constraints they can place on the ag es, distances, and 
progenitor mas ses of magnetars (e.g.. IVrba et al.l 120001 
IGaensler et "all 120011; iMuno et al.ll2006f) We discuss this 
in more detail below. 

The structure of this paper is as follows. First, in § 11.11 
we discuss the two targets of this analysis, SGR 1900+14 
and IE 2259+586. We then discuss the archival data 
and our new Chandra X-ray Ob servatory Advanced C CD 
Imaging Spectrometer (ACIS: iGarmire et al.ll2003T) ob- 
servations in § [2] In § [31 we present a detailed analysis of 
those data, along with archival data used to test our tech- 
niques and simulations. We present our proper motion 
measurements in §01 along with comparison to other sim- 
ilar measurements and the prospects for improvement. 
Finally, we give our discussion and conclusions in § [5j 

1.1. SGR 1900+14 and IE 2259+586 

The third S GR discovered 

(|Mazets. Golenetskii. fc Gurvanl [l979h . SGR 1900+14 
was identifi e d as a persistent X-ray source by 
iHurlev et alJ (|1996l ). Subse quent observations iden- 
tified the 5.2-s pulse pe riod (IHurlev et~aT. [l999c) with 
significant spin-down (jKouveliotou et aL 1 199£ ) im- 
plying a large dipolar magnetic field. The magnetar 
nature of this source was confirm ed by the giant flare 
detected on 1998 August 27 ([Hurley et alJ Il999at 
iThompson fc Duncadll995l. Il996h. 

Vasisht et al.l(|1994D and lKouveliotou et all (|1994l ) sug- 
gested that SGR 1900+14 could be associated with 
the nearby supernova remnant (SNR) G42.8+0.6, based 
largely on position coincidence and the possible associa- 
tion of a number of other AXPs and SGRs with SNRs. 
However, lLorimer fc Xilourii ([2000) called the associ- 
ation into question when they discovered a relatively 
young (10 4 -10 5 y r old) radio pul s ar nea r both the SGR 
and the SNR, and IGaensler et all (|2001f ) calculated a 4% 
chance of random alignment between the SGR and SNR. 
This chance is even greater if one accepts that super- 
nova explosions would naturally tend to be clustered, 
so SGR 190 0+14 appearing nea r SNR G42.8+0.6 (and 
other SNRs; Kaplan et all 12002) is not unlikely. If the 
association between SGR 1900+14 and SNR G42. 8+0.6 



were real, SGR 1900+14 would need a significant trans- 
verse velocity to make it from the explosion site to its cur- 
rent lo cation outside the SNR's radio shell ([Hurley et al.l 
I1999H) . SGR 1900+14 is w 15' from the center of the 
SNR (at approximately a J20 oo = 19 h 07 m 01 s , <5j 2 ooo = 
+09°04'10" from lo w- frequency radio i mages, with an 
uncertainty of ±30"; [Kaplan et"al"1l2002f ). This implies a 
proper motion of 90t 1 ~q L masyr" 1 or a transverse velocity 
v± « 2100d5T 10 1 kms -1 , where the distance to the sys- 
tem is 5d5 kpc, and the age 6 is I Otiq kyr. This is a large 



velocity for a neutron star 



not unheard of ([Chatteriee et al 



Hob bs et al.l 120051 ). but 



2005h . 



An alternati ve associat i on ha s been proposed for 
SGR 1900+14. IVrba et al.l (12001 identified a cluster of 
massive stars, with several M supergiant members and 
> 10 other stars, about 12" away from the SGR 1900+14. 
This has been viewed as a much more secure association 
for SGR 1900+14, since suc h clusters are r a rer an d it is 
much closer on the sky, and Wachtcr ct al. (2008) claim 



5 IHelfand et al.l 1120071 ) also used infrared measurements of 
XTE J1810— 197 to constrain the proper motion, and found a much 
coarser but consistent result. 



that the infrared ring found around SGR 1900+14 could 
only be powered by the cluster stars. The possible in- 
consistency betw^ejitlieejrtmction to SGR 1900+14 and 
the cluster (Kaplan ct al. 2002) may no t be significant, as 
demonstrated bv lWachter et alJ ([2008D . This cluster is at 
12-15 kpc, so if the association is real then SGR 1900+14 
is at a distance of 4 ~ 3. 

The anomalous X-ray pulsar IE 2259+586 was the 
first member of its cla ss to be discovered, when 
IFahlman fc Gregory! (|1981f ) identified it as an X-ray point 
source in the center of the SNR CTB 109 that had co- 
herent 7-s pul sations. The associa tion with the SNR 
seems secure (Gacnsle r et alJ [20011 ). as the positional 
coincidence is very good (< 3') compared to the lo- 
cal density of SNRs. Based on interactions between 
the SNR and H II regi o ns wi th measured distances, 
iKothes. Uvaniker. fc Yarl ([2002) find a distance 7 of « 
3 kpc. 

2. OBSERVATIONS 

We examined the Chandra archive for suitable obser- 
vations of magnetars. We found data on SGR 1900+14 
and IE 2259+586 that we used for the first epochs of 
our proper motions studies. These data were taken with 
the ACIS instrument in the full-frame mode and with no 
grating, in contrast to many other observations of magne- 
tars that use subarrays, continuous clocking, or gratings 
to get better timing or spectral resolution. 

We re-observed SGR 1900+14 and IE 2259+586 with 
Chandra, with the new data serving as second epoch ob- 
servations to try to measure or constrain their proper 
motions. For the second epochs, we chose to gener- 
ally replicate the first epochs in detector choice, aim- 
point, and exposure time. One of the priorities was 

6 Both the age and distance of the SGR 1900+14 are highly 
uncertain. The spin-down age P/2P is about 1 kyr, but the spin- 
down rates of magnetars can easily vary by factors of several 
( Woods & Thompson 2006), thus making the spin-down age an un- 
reliable estimator of the true age. The distance to SGR 1900+14 
itself is largely unconstrained. 

7 There is some uncertai nty about the association and he nce the 
distance to IE 2259+586. Dur ant fc van Kerkwiikl {2006) find a 
distance to IE 2259+586 of 7.5 ± 1.0 kpc using the "red clump" 
method. They argue that CTB 109 could also plausibly be at this 
distance too, based on a reinterpretation of radial velocity data, 
but there are uncertainties associated with this analysis. 
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TABLE 1 

Summary of Chandra Observations 



Target 


Epoch 


ObsID 


Date 


MJD 


Aimpoint 


Exposure 


Rotation 














(ksec) 


(deg) 


SGR 1900+14 


1 


1954 


2001-Jun-17 


52077 


ACIS-I 


29.8 


148.6 




2 


6731 


2006-Jun-04 


53891 


ACIS-I 


24.6 


134.0 


IE 2259+586 


1 


725 


2000-Jan-12 


51555 


ACIS-S 


18.9 


312.8 




2 


6730 


2006-May-09 


53864 


ACIS-S 


24.8 


73.7 



Note. — The rotation given is the nominal roll angle of the spacecraft, defined by the 
R0LLJJ0M header keyword. 



to detect as many other sources as possible besides the 
SGR 1900+14/1E 2259+586: as discussed below, our 
accuracy is limited by the number and quality of the ref- 
erence sources. In Table[T]we give the details of both the 
archival observations and our new ones. Note that the ro- 
tations of the new and old observations are in general dif- 
ferent. For SGR 1900+14, which used the square ACIS-I 
detector, this made no difference, but for IE 2259+586 
the reference sources changed somewhat between epochs, 
giving us a smaller set than we would have liked: see dis- 
cussion in § 14.21 

3. X-RAY ASTROMETRY 

Our general analysis technique is to (1) measure the 
positions of a number of field X-ray sources ("refer- 
ence sources") common to both epochs; (2) measure 
the position of the object of interest (SGR 1900+14 
or IE 2259+586) in both epochs; (3) use the refer- 
ence sources to determine a transformation between the 
two observations; and (4) apply that transformation to 
the object of interest, finding the position difference (if 
any) between the two epochs that we can then convert 
to a pro per motion. This g e nerally echoes the tech - 
niqu es oflHui fc Becker! (120061) , iWinkler k. Petrel (120071) . 
and IMotch et all (|2007al lb[). although like iMotch et al l 
(|2007al lbD we prefer to use many reference sources and 
not just a few. 

We must determine the best way to carry out each 
step, and assess the uncertainties associated with each. 
For this reason, we discuss each step at length below, in 
many cases using data from other observations and/or 
simulations to assess accuracy. 

In the discussion below, we work as much as possi- 
ble in the sky pixel (£ s ky,2/sky) frame. This is a frame 
constructed by the Chandra event processing that ro- 
tates the detector (xdet, Diet) frame to have y s \ y point- 
ing north and a; s k y pointing west. While in general we 
prefer to stay in a frame with as little processing as possi- 
ble, which would argue for retaining the detector frame, 
the continuous dither of Chandra (in a Lissajous pattern 
with amplitude ±8" and periods of 707 s and 1000 s along 
the two axes) makes that difficult. The conversion from 
detector to sky frame involves correcting for that dither, 
for the nominal orientation of the spacecraft, and for any 
other position excursions, but this aspect reconstruction 
is of sufficient quality that it does not limit our accuracy. 

3.1. Reference Sources 

For our main position-fi nding algorithm, w e chose 
the CIAO task wavdetect (Freeman et al. 2002), which 
uses a wavelet-based detection algorithm. This echoes 
the choices made by such large X-ray surveys as the 



Chandra Multiwa velength Plane (ChaMPlane) Survey 
(Hon g~et alj 120051) and the Cha ndra Multiwavelength 
Project (ChaMP: lKim et al.ll2007l) . In those cases the al- 
gorithm choice was made for optimum source detection, 
with astrometry only as a secondary goal, but they did 
find that wavdetect produced the most reliable astrom- 
etry. We verified this choice by comparing the results 
of wavdetect, the sliding-cell task celldetect, and our 
own implementation of an interactive centroiding task 
that operates on the raw event lists, and found that 
wavdetect gave the smallest scatter in the various tests 
described below. 

As we were not overly concerned with detecting every 
possible source or measuring such source properties as 
count-rate, we ran wavdetect on the raw level 2 event 
files returned from the Chandra archive without any ad- 
ditional processing (e.g., exposure map creation), using 
wavelet scales of 1,2,4,8 pixels. We did test whether 
removing sub-pixel randomization 8 made a difference in 
the positions, and we find no systematic difference. 

However, the uncertainties reported by wavdetect are 
systematically low, especially for off-axis sources. We 
therefore investigated the uncertainty models used by 
the ChaMP a nd ChaMPlane pr oject s. Independent o f 
the analyses of iHong et ail (|2005h and lKim et ail (|2007h . 
which relied on MARX (Model of AXAF Response to X- 
rays) and SAOsac (a raytracing code developed for Chan- 
dra) simulations of sources, we wished to see how reliable 
the position measurement and uncertainty estimation of 
wavdetect are. Therefore we took a number of relatively 
long ACIS-I exposures from the archive with a large num- 
ber of point sources. We then split those exposures into 
two halves by time, effectively making two sub-exposures 
that are otherwise identical, and processed each with 
wavdetect in the same manner. This way we can exam- 
ine how the positions from the first half agree with those 
from the second half: since the two sub-exposures are 
from the same observation, the positions should agree. 
Of course the sources in each sub-exposure with have 
fewer counts than the sources in the total exposure, but 
there are still a range of count levels and off-axis angles. 

We first did this comparing the results of wavdetect 
with celldetect and an interative centroiding algorithm 
that we implemented, and found wavdetect to have the 
lowest dispersion. We then varied the energy range of our 
source extraction and found that including event energies 
from 0.3 keV to 7 keV optimized the signal-to-noise of 
the extraction. 

We then compared the position differences between 
the first and second sub-exposures with the position 

8 See http://cxc.harvard.edu/ciao/threads/acispixrand/ 
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uncertainties returned by wavdetect. We found that 
wavdetect significantly underestimated the uncertain- 
ties, especially at off-axis angles 9 > 3.5'. The com- 
parison gave a reduced x 2 of 4.31 for 1450 degrees-of- 
freedom (DOF), and the shape of the distribution is 
clearly broader than the expected distribution, with more 
points at > 2a. In contrast, both the ChaMP and ChaM- 
Plane models work well, giving reduced \ 2 values of 1.33 
and 1.02. The uncertainty from the ChaMPlane model 
exceeds that from the ChaMP model slightly at large and 
small off-axis angles, leading to the slightly lower \ 2 f° r 
ChaMPlane. 

Overall, both the ChaMP and ChaMPlane uncertainty 
models appear to work. Both are based on simulations 
using MARX, and so the agreement is not surprising, but 
since they agree with each other and also agree with the 
actual data presented here (see § 13. 4[) we feel comfortable 
using either. The functional form of the ChaMP model is 
somewhat simpler and conforms more to the expected re- 
lation between uncertainty a XtV and signal-to-noise ratio, 
°- 46 (C is the number of counts above the 



with a x , y oc C 



background) compared to the expected a XA 



C 



-0.5 



(i.e. 



a x>v oc (S/N)^ 1 / 2 ). In contrast, the ChaMPlane model 
is a more complicated function of C involving powers of 
the logarithm of C. Therefore we w ill use the ChaMP 
model (see eq. 14 of lKim et aLll2007[ ): 



a X;V = 0.77 CT °- 46 10 one pixels (1) 

with a x y as the x or y uncertainty (con verted from 
the radial un certainty in [Kim et ah |2007| by di viding 
by 1.515; see lLampton. Margon. fc Bowverl fl976[ K and 
6 in arcminutes. Note that this is only valid for low 
count levels (C < 131) as is the case here but not for 
SCR 1900+14 or IE 2259+586. 

3.2. SGR 1900+14 

For SGR 1900+14, we are in a very different regime 
than for the reference sources. Instead of relatively faint 
sources at a variety of off-axis angles, we have a bright 
source at 6 < 0.5'. A different approach to position mea- 
surement and uncertainties than that described in § 13.11 
is needed. 

The ACIS-I images of SGR 1900+14 have the source 
at an off-axis angle of 20", and with about 5000 counts. 
Technically, then, the ChaMP uncertainty models do not 
apply (they are only valid for < 2000 counts). Even so, 
the ChaMP model for high numbers of counts does not 
agree with what we would naively expect, with a x ^ y oc 
C~ 02 . We therefore did our own SAOsac/MARX (using 
the Chandra Ray Tracer ChaRT as an inte rface to the 
SAOsac raytracing code; iCarter et al.ll2003r ) simulations 
of this source. 

We simulated a bright point-source at the cor- 
rect off-axis location with the approximate spectrum 
of SGR 1900+14: a power-law with photon index 
r = 1.9 and absorpti on with TVh = 2 x 10 22 cm~ 2 
dMeregfietti et all 12006ft . We simulated a source with 
700,000 counts, many more counts than are actually de- 
tected in our exposures of SGR 1900+14, so that we 
could divide this exposure into sub-exposures as above 
for comparison. We made a series of sub-exposures, di- 
viding the 700,000 total counts into 12 sub-exposures 
with 56,000 counts, 25 sub-exposures with 27,000 counts, 




0.5 1 1.5 2 2.5 3 

Radius (pixels) 

Fig. 1. — Radial profile of SGR 1900+14 (points) compared to 
that of a MARX simulation with 30 times the number of counts (solid 
line) that does not include the effects of pileup. The agreement 
between the data and the model illustrates the small effect of pileup 
on the radial profile of SGR 1900+14. 

and so on, down to 300 sub-exposures with 2200 counts. 
We then ran wavdetect on each of the sub-exposures in 
a series and compared all of the positions. In contrast 
to the situation with the reference sources, where the 
raw wavdetect uncertainties underpredicted the uncer- 
tainties, here they did well. In all cases the comparisons 
yielded reduced \ 2 values consistent with 1.0. Fitting the 
uncertainty as a function of number of counts, we find a 
relation close to the expected one, with a x>y oc C~ 0A7 . 

To confirm this, we also divided up the real exposures 
of SGR 1900+14 into 10 sub-exposures of sa 500 counts 
each. Each position as measured by wavdetect had an 
uncertainty of 0.04 pixels, and again the positions of the 
sub-exposures were all consistent with each other within 
the uncertainties. 

A final concern is pileup. With a count-rate of ~ 
0.17s -1 , the expected pileup fraction is ~ 20%, with 
~ 0.5 event expected in each 3.2 s frame. This means 
that our data are affected by pileup such that spectral 
estimation will not be robust, but at this pileup level the 
effect on astrometry is minimal. The spatial distribu- 
tion of events is still largely consistent with the point- 
spread function (Fig.[l|. Our MARX simulations here did 
not account for pileup, but even so they were consistent 
with our examination of the sub-exposures of the real 
data. Overall, then, wavdetect positions and uncertain- 
ties seem sufficient for observations of SGR 1900+14. 

3.3. IE 2259+586 

For IE 2259+586, we could not use the results of 
wavdetect. The pileup for this source is severe enough 
that the X-ray image has a central hole, where so many 
X-ray photons landed that all were judged to be cosmic 
rays. We also have issues at large radii: beyond a radius 
of ks 4", there is extended emission from a combination 
of the SNR in which IE 2259+ 586 is embedded as well 
as a halo from dust scattering (jPatel et al.ll200lT ). 

We initially tried doing simple centroids on the event 
data, where we used events in an annulus whose inner ra- 
dius varied from to 5 pix (to eliminate the central hole) 
and whose outer radius varied from 15 pix to 100 pix. For 
data with no central hole, iterating this process will con- 
verge on the correct center (this was one of the methods 
we used in 8 13.21). However, we found that this method 
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TABLE 2 
Source Positions for SGR 1900+14 



ID Epoch 1 Epoch 2 





•''sky 


11 l 

.ysky 


CJounts 




■^sky 




ysky 




(Jounts 


^circiriinj 


X 


4075.5 ±0.0 


4067.4 ± 0.0 


5232.3 


0.3 


4109.7 ±0 





4022.9 ±0 





3805.6 


0.3 


1 


3991.4 ±0.5 


4425.6 ±0.5 


15.5 


2.8 


4025.8 ± 


5 


4380.5 ± 


5 


11.7 


2.8 


2 


4080.5 ±0.4 


4494.6 ± 0.4 


30.5 


3.3 


4114.7 ±0 


6 


4448.8 ± 


6 


12.7 


3.2 


3 


4159.6 ±1.0 


4707.2 ± 1.0 


10.5 


5.0 


4192.8 ± 1 


1 


4662.8 ± 1 


1 


8.1 


5.0 


4 


4385.9 ±0.4 


4822.8 ±0.4 


138.2 


6.4 


4419.6 ±0 


5 


4779.3 ± 


5 


95.3 


6.4 


5 


4409.6 ±0.7 


4421.9 ±0.7 


9.6 


3.7 


4442.6 ± 1 


3 


4378.2 ± 1 


3 


2.9 


3.7 


6 


4025.4 ±1.5 


4785.1 ± 1.5 


6.3 


5.7 


4059.5 ± 1 





4743.5 ± 1 





15.8 


5.7 


7 


4383.0 ±1.7 


4655.8 ± 1.7 


3.5 


5.2 


4418.8 ± 1 


1 


4612.8 ± 1 


1 


5.7 


5.2 


8 


4395.0 ± 1.9 


5093.7 ± 1.9 


18.6 


8.5 


4425.4 ± 2 


3 


5053.4 ± 2 


3 


12.5 


8.6 


9 


3638.7 ±0.6 


4492.5 ± 0.6 


29.0 


5.0 


3674.1 ±0 


8 


4446.7 ± 


8 


15.8 


4.9 


10 


3997.8 ±0.5 


4280.2 ± 0.5 


5.7 


1.7 


4032.3 ±0 


5 


4235.4 ±0 


5 


6.8 


1.7 


11 


3055.8 ±1.9 


4475.3 ± 1.9 


23.9 


9.1 


3090.0 ± 1 


9 


4426.7 ± 1 


9 


24.1 


9.0 


12 


3112.3 ± 1.0 


4680.9 ± 1.0 


128.4 


9.4 


3148.0 ±2 





4635.6 ± 2 





24.8 


9.3 


13 


3319.6 ±1.3 


4314.9 ± 1.3 


13.8 


6.6 


3355.5 ± 1 


2 


4270.0 ± 1 


2 


16.7 


6.5 


14 


4313.9 ±0.4 


4043.5 ± 0.4 


9.7 


1.8 


4348.5 ±0 


5 


3999.2 ± 


5 


7.7 


1.9 


15 


4398.5 ±0.6 


4043.0 ± 0.6 


7.7 


2.5 


4433.8 ± 


7 


3998.2 ± 


7 


5.8 


2.6 


16 


4706.4 ± 0.9 


3894.6 ± 0.9 


14.0 


5.3 


4739.4 ± 1 


1 


3848.5 ± 1 


1 


8.9 


5.3 


17 


4781.0 ±0.7 


3963.5 ± 0.7 


35.0 


5.7 


4814.3 ± 1 





3918.2 ± 1 





15.3 


5.8 


18 


4694.9 ± 1.1 


4248.4 ± 1.1 


9.1 


5.1 


4730.5 ± 1 





4200.2 ± 1 





10.1 


5.1 


19 


3839.5 ± 0.6 


3917.6 ±0.6 


7.7 


2.6 


3873.1 ±0 


5 


3872.6 ± 


5 


10.9 


2.5 


20 


3931.1 ±0.5 


3957.6 ±0.5 


8.7 


1.8 


3965.4 ±0 


1 


3912.9 ±0 


1 


9.8 


1.7 


21 


4005.7 ±0.4 


3945.3 ±0.4 


8.6 


1.4 


4040.1 ±0 


•1 


3900.6 ± 


1 


8.7 


1.4 


22 


4046.3 ±0.3 


3940.9 ±0.3 


13.6 


1.3 


4080.3 ± 


5 


3896.7 ± 


5 


5.8 


1.3 


23 


3582.5 ±0.9 


3748.2 ± 0.9 


12.9 


5.1 


3617.3 ± 1 





3703.3 ± 1 





10.4 


5.0 


24 


3925.2 ± 1.1 


3545.0 ± 1.1 


7.4 


4.7 


3958.9 ±0 


8 


3499.8 ± 


8 


13.9 


1.7 


25 


4125.5 ± 1.0 


3069.3 ± 1.0 


72.4 


8,1 


4158.6 ± 1 


1 


3024.0 ± 1 


1 


53.8 


8.5 


26 


4207.8 ±1.4 


3368.5 ± 1.4 


8.5 


6.0 


4241.9 ± 1 





3323.4 ± 1 





16.3 


6.1 


27 


4255.9 ±0.9 


3715.3 ±0.9 


4.8 


3.1 


4289.4 ± 


7 


3670.6 ± 


7 


7.8 


3.1 



Note. — Positions are in sky pixel coordinates, 6 is the off-axis angle, and source X is SGR 1900+14. 
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Fig. 2. — Radial profiles of IE 2259+586 (squares and circles 
for epochs 1 and 2) compared to that of a MARX simulation that 
includes pileup (x's). We fit the data outside of a radius of 4pix 
(vertical dotted line) with a Gaussian (dashed line) multiplied by 
a hole (dot-dashed line; the hole has the form tanh(r/ro) 10 with 
ro = 1.25 pix); the full model is the solid line. 

gave results that did not converge: changing the outer ra- 
dius led to systematic changes in the centroid position. 
This was partially due to asymmetries in the extended X- 
ray emission, but more important were the effects of non- 
uniform exposure over the AC1S CCDs. Together, these 
effects led to variations (both smooth and discrete) of the 
overall count level over the image, and when the annuli 
that we chose included those portions, they skewed the 
resulting centroid. Simulating the data in MARX, where 
we used similar exposure maps, reproduced the wander- 
ing centroid, although not to a high enough precision that 
the MARX data would all ow us to correct the obs ervations. 
Instead, we followed iHulleman et al.l (|2001l ) and ex- 



plicitly used a local, symmetric function to find the po- 
sition of IE 2259+586. This is similar to what we did 
with SGR 1900+14, where wavdetect found the posi- 
tion that maximized the overlap between the data and 
a wavelet function, but for IE 2259+586 we use our 
own kernel function and our own routines to find the 
position. For the kernel we used a Gaussian function 
multiplied by a hole, meant to mimic the shape of the 
piled-up image. The function was radially symmetric. 
IHulleman et al.l (|2001[ ) used a hyperbolic tangent for 
their "hole" function: we found that this did not de- 
cline quite quickly enough to match our data, and in- 
stead used tanh(r/ro) , where ro is the radius scale and 
the exponent serves to increase the sharpness of the hole; 
the results are not sensitive to the exponent as long as 
it is > 4. A fit to the data from both epochs (averaged 
in azimuth), along with a MARX simulation that includes 
pileup, is shown in Figure [2 To make that fit, we fit the 
Gaussian to the data, the MARX model with pileup, and 
the MARX model without pileup, all for radii of > 4 pix. 
At those radii the data and the two models all agree rea- 
sonably well. Within 4 pix the effects of pileup become 
more severe, and the MARX model is not able to accurately 
predict the data, but the agreement is qualitatively rea- 
sonable. We found a hole radius r « 1.25 pix, but again 
the fit was not extremely sensitive to this value. The goal 
of this model is to give a good qualitative representation 
of the data, not to fit it in great detail. 

To actually find the position of IE 2259+586, instead 
of fitting the data, which is sensitive to the details of the 
model as well, we prefer to cross-correlate a fixed model 
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TABLE 3 
Source Positions for IE 2259+586 



ID 




Epoch 1 








Epoch 2 






•^sky 


Vsky 


Counts 


(arcmin) 


■^'sky 


J/sky 


Counts 


(arcmin) 


Y 


4141.7 ±0.1 


4171.9 ±0.1 




0.7 


4117.0 ±0.1 


4070.7 ± 0.1 




0.3 


1 


4093.8 ±0.4 


3567.5 ±0.4 


64.6 


4.3 


4069.9 ±0.3 


3466.8 ± 0.3 


157.7 


5.2 


2 


4536.1 ±0.2 


3899.2 ±0.2 


295.7 


4.0 


4510.9 ±0.5 


3798.9 ±0.5 


26.6 


4.2 


3 


4159.5 ±0.3 


3799.6 ±0.3 


38.7 


2.5 


4135.0 ±0.4 


3699.7 ±0.4 


21.6 


3.3 


4 


4270.1 ±0.2 


4755.3 ±0.2 


330.9 


5.6 


4245.3 ±0.3 


4654.4 ±0.3 


103.4 


1.7 


5 


4350.6 ±0.4 


4446.1 ±0.4 


28.9 


3.5 


4326.3 ±0.3 


4344.7 ±0.3 


47.6 


2.8 


6 


4273.3 ± 1.2 


4961.5 ± 1.2 


24.4 


7.2 


4246.7 ± 1.0 


4857.9 ± 1.0 


20.8 


6.1 



Note. — Positions are in sky pixel coordinates, 9 is the off-axis angle, and source Y is IE 2259±586. 
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Fig. 3. — Cross-correlation of the first epoch of data for 
IE 2259±586 (black solid contours) with the model from Eqn. [2] 
(red dashed contours). The model and data have been shifted to 
the origin based on the result of the cross-correlation. The results 
for epoch 2 are similar. 

with the data. Our fixed model is that shown in Figured) 



Fir) = e- r2/2a ' 2 tanh 
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'o 



(2) 



where we found a = 2.5 pixels and tq = 1.25 pixels gave 
good results. This function is radially symmetric, and by 
cross-correlating rather than fitting we can ignore issues 
of normalization and the variations in the wings of the 
radial profile. Above, we could have just fit the data to 
a Gaussian function only for radii greater than 4 pixels, 
but this would not have allowed the cross-correlation. 
Using a Gaussian kernel function with no hole but just 
excluding a portion of the data performed poorly, since 
the sharp "horns" of the data outside of the hole can 
cause local maxima in the cross correlation if the kernel 
function has a single peak like a Gaussian. 

We cross-correlated images of the filtered event files 
created from both the data and the MARX simulations 
with the model (Fig. [3]). The cross-correlation results 
were smooth and well-behaved, and we interpolated be- 
tween our grid points to find the best-fit offset. We made 
sure that the fitted position did not depend on the de- 
tails of the cross-correlation: we varied the binning, the 
outer radius, and the parameters of the model (er and r ), 
and in all cases the resulting variations in the position 
were < 0.05 pix. As a further check, we shifted the event 



positions by random amounts of up to 1 pix, and found 
that we were able to recover the shifts to an accuracy 
of < 0.1 pix; this test also verified that our result did 
not depend on the origin of the binning. Based on these 
tests, as well as additional MARX simulations, we estimate 
that our cross-correlation is accurate to ±0.1 pix in each 
coordinate. 

As another test, we also cross-correlated the binned im- 
ages from epochs 1 and 2 against each other, rather than 
against a smooth kernel. In general the results of this test 
were similar, where we could recover an input shift at a 
level of < 0.1 pix, and the offsets agreed to a similar pre- 
cision. We plot the results of all cross-correlations (epoch 
1 vs. model, epoch 2 vs. model, and epoch 1 vs. epoch 
2) in Figure [H where we plot the results relative to the 
best-fit positions given in Table [3] The cross-correlations 
that we plot do not agree exactly with the positions in 
Table |3] as those positions were derived using multiple 
cross-correlations with random offsets used to eliminate 
the effects of binning, but the agreement is better than 
the quoted ±0.1 pix uncertainty. 

For a final check, we examined the position of the read- 
out streak on the two epochs. This streak comes about 
because the ACIS detector has no shutter, so the CCDs 
receive photons even as the charge is being clocked off. 
This gives rise to a one-dimensional streak centered on 
each bright source which is aligned with the CCD row 
and column axes. It ends up having ~ 2% of the to- 
tal number of events as the original source 9 , and can be 
used for both positional and spectral information. Here 
we make use of only the former. 

To use the streak, we took the event positions from 
the level-2 event file. We then rotated by the R0LLJJ0M 
header keyword, putting them into a frame where the 
CCD rows and columns were vertical and horizontal 
(similar to the CHIP frame). The readout streak, which 
is a line in the sky (x,y) frame not aligned with either 
axis, is now horizontal. We then stepped along the line, 
computing the median event position in the vertical di- 
rection for bins of 50 pix width in the horizontal direc- 
tion: see Figure [5j For both epochs 1 and 2, we found 
streak positions that were consistent to ±0.02 pix with 
our best-fit positions of IE 2259±586 from Table[3j albeit 
in one dimension. While not sufficient to get a position 
for IE 2259±586 alone, this provides additional confir- 
mation that our cross-correlation technique is accurate 
to the claimed 0.1 pix. 

9 See http : //exc .harvard. edu/cal/Acis/Cal_prods/xfer_streak/index. html 
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Ep. 1 vs. Ep. 2 

Ep. 1 vs. Ep. 2 (y) 

Ep. 1 vs. Model {x}-0.05 

Ep. 1 vs. Model (y}-0.05 

- - Ep. 2 vs. Model (x)-0.1 

- - Ep. 2 vs. Model (y)-0.1 




-0.5 0.5 
Position Offset (pixels) 

Fig. 4. — Cross-correlation results for IE 2259+586. We show 
one-dimensional cuts through the peak of the cross-correlation sur- 
face, which we use to identify the best-fit positions in x (blue 
curves) and y (red curves). For each cross-correlation (epoch 1 vs. 
epoch 2: dot-dashed line; epoch 1 vs. the model: solid line; epoch 
2 vs. the model: dashed line; in all cases the model is from Eqn.[2j 
we interpolated to find the maximum value of the cross-correlation, 
which we identify as the best-fit position, and the cross-correlations 
have been offset vertically for clarity. The data plotted here have 
had our nominal positions subtracted, giving offsets near 0. They 
are not quite at zero as the nominal positions were the averages 
of cross-correlations done with a range of random offsets added to 
mitigate the effects of binning, but the differences are consistent 
with our claimed uncertainty of 0.1 pix (vertical dotted lines). 



Displacement along streak (pixels) 
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-300 -100 100 300 500 

Displacement along streak (pixels) 

Fig. 5.— Fits to the readout streak of IE 2259+586: epoch 
1 (top) and epoch 2 (bottom). We rotated the event locations 
according to the R0LLJJ0M header keyword, ending up in a frame 
where the CCD readout streak was horizontal. The grayscales are 
the binned images for each epoch; the bins appear rectangular in 
this plot because of the drastically different scales for the ordinates 
and abscissae. We then constructed 50 pixel wide bins along the 
streak and computed the median event locations perpendicular to 
the streak, which we show as the points. The dashed lines are the 
averages of the median positions. We have excluded events within 
a radius of 20 pix around IE 2259+586. 

3.4. Frame Transformations 

We considered a number of ways to transform between 
the position measurements at each stage. To test this, 
we used 13 ACIS-I observations 10 of the Hubb le Deep 
Field North fHDF-N; lHornschemeier et al.ll2000l and sub- 
sequent papers). The observations were between 50 ks 
and 140 ks in depth, and detected > 100 sources. Note 

10 We used ObsIDs: 2232, 2233, 2234, 2421, 2423, 3293, 3294, 
3388, 3389, 3390, 3391, 3408, and 3409. 



that the observational setup here was, if anything, more 
taxing that what we used for our proper motion obser- 
vations. This was because a number of the observations 
had a different aimpoint (off by about half of an ACIS de- 
tector) and varying roll angles. We identified 38 sources 
that were common to each exposure (note that this was 
not a completely exhaustive search, but comprised those 
that we could identify easily) and computed transfor- 
mations between each pair of observations. We used a 
number of different transformations with different levels 
of complexity. 

2- parameter : A shift in x and y. 

3- parameter : A shift in x and y, along with a rotation. 

4- parameter : A shift in x and y, along with a rotation 

and an overall scaling. 

6-parameter : A general linear transformation, in- 
volving a shift and a separate rotation and 
scaling for each axis (e.g., see Eqn . A5 in 
IRaplan, van Kerkwiik. & Anderson|[2007l ). 

We were transforming between x and y values in the 
SKY frame, so the nominal position angle of Chandra 
was already taken out. Any remaining rotation or scal- 
ing would account for imprecise knowledge of the space- 
craft's orientation and breathing of the ACIS-I detector 
(changes in the plate-scale due to thermal variations), 
respectively. 

Overall we found that the 3-parameter transforma- 
tion provided a significant decrease in \ 2 compared to 
the 2-parameter transformation, with a mean A% 2 of 
23.8 going from 74 to 73 DOF. This is extremely sig- 
nificant, with an F-test probability of < 10 -5 that this 
is due to chance. Therefore a non-zero rotation is re- 
quired between different observations, with a magnitude 
of « 0.05°. Comparing sets of three observations, the ro- 
tations cf> fitted between all of them are consistent, with 
(e.g.,) 0(1 -> 2) + 0(2 ->• 3) » 0(1 -> 3), indicating that 
the rotations can be physically meaningful. With just the 
3-parameter transformation the reduced % 2 was consis- 
tent with 1.0, and the distribution of position residuals 
was consistent with the expected distribution (Fig. [5]). 
This indicates that our uncertainty model is reasonable: 
tests using the ChaMPlane model gave similar results, 
while those using the raw wavdetect uncertainties had 
reduced \ 2 larger than 2. Beyond this, the more compli- 
cated transformations were not warranted, with a median 
A% 2 of only 4.2 between 3 parameters and 6. 

4. MEASURING PROPER MOTIONS 

4.1. SGR 1900+14 

To actually measure the proper motion of 
SGR 1900+14, first we ran wavdetect on both ob- 
servations. The parameters were the same as those 
discussed in § [3j We then identified sources common to 
both epochs of the data, finding 27 sources in addition 
to SGR 1900+14. These sources were on all four of 
the ACIS-I CCDs. It would be preferable to have the 
SGR 1900+14 centered on a single CCD, and to only use 
sources on that detector. That way we would not have 
to worry about inaccuracies in the relative positions 
of the detectors or shifts in their positions with time 
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Fig. 6. — Position differences for 2964 sources compared between 
different observations of the HDF-N (38 unique sources and 78 pairs 
of observations), normalized by position uncertainty (computed ac- 
cording to Eqn.[TJ- The data are the histogram, and we also show 
the expected distribution with N(r) oc re~ r I' 2, for comparison. 

affecting our astrometry. However, that would severely 
restrict our field of view, and we felt that the gain of 
additional reference sources outweighed the possibility 
of smaller uncertainties for the detected sources. The 
various tests and discussion in § 13.11 seem to indicate 
that using sources on different CCDs does not adversely 
affect our astrometry, at least at the level considered 
here. 

With the source detections, we then updated the 
position uncertainties of all of the sources except 
SCR 1900+14 using Eqn.[Q Finally, we fit for the trans- 
formation between the two epochs of data, using a 3- 
parameter fit (§ 13. 4| . The data gave a very good fit, 
with x 2 — 37.7 for 51 DOF (27 reference sources and 
3 parameters). There were no significant outliers in the 
data, which is not really surprising since we only included 
sources that we could match initially between the epochs. 
The magnitude of the residuals is ~ 100 masyr -1 , which 
is far larger than the proper motion one would expect 
from random field stars (i.e., at 1 kpc lOOmasyr -1 is 
475 km s , which is much larger than the velocity dis- 
persion of Galactic stars). Therefore we can consider that 
we are measuring the motion of SGR 1900+14 relative 
to reference stars that are at rest. Similarly, corrections 
for Galactic rotation and Solar motion are too small to 
be considered. 

The fitted rotation here was consistent with 0. We 
therefore found that a 2-parameter transformation was 
sufficient (with virtually no change in x 2 ), and with the 
smaller number of free parameters the resulting position 
difference had smaller uncertainties. While in § 13.41 we 
found that a rotation was necessary, here we do not, 
although whether it is the simpler observational setup 
(both observations had the same aimpoint, in contrast 
to § I3.4p or just the reduced accuracy from a smaller 
number of reference sources that changed the situation 
we do not know. We plot the residual position differences 
of the reference sources in Figure [7] 

As a check on our transformation, we performed jack- 
knife tests, where we calculate the transformation for 
the data excluding a single reference source (see, e.g., 
lEfron fc Tibshiranil [l993h . This tests whether or not 
the transformation depends overly on a single reference 



source, and indeed we find that it does not. We find 
that the means of the jackknifed shifts are consistent 
with the value using all of the data, and that the stan- 
dard deviations are significantly less than our estimated 
shift uncertainties. We then performed a bootstrap test. 
Here we created 100 data sets, each including a random 
subset chosen with replacement from the 27 reference 
sources. The resulting shifts are again distributed around 
the proper values created from the actual data, and the 
standard deviations are comparable to our uncertainties 
(they are actually slightly smaller). Based on these two 
tests, we see that our transformation calculation is rea- 
sonably robust to choices of the input data. 

With the 3-parameter transformation, we find that 
the position of SGR 1900+14 is consistent between the 
two epochs, moving by 0.27 ± 0.27 pixels, implying 
no measurable proper motion. Note that the bulk of 
this uncertainty comes from uncertainties in the fitted 
transformation (uncertainties in each coordinate of the 
shift of 0.26 pixels, and an uncertainty in the rotation 
of 0.03°), since the individual position measurements 
of SGR 1900+14 had uncertainties of < 0.02 pixels. 
We also checked to see if the choice of transformation 
(2, 3, 4, or 6 parameters), uncertainty model (ChaMP 
or ChaMPlane), or level of processing (with or with- 
out ACIS randomization) made a difference, and find 
that all of the results were consistent. With the 2- 
paramcter transformation, which as discussed above is 
sufficient, we find a very similar shift but smaller resid- 
uals: (Ax, Ay) = (0.02 + 0.19,0.26 + 0.19) pixels, or 
a magnitude 0.27 ± 0.19 pixels for the 2-parameter fit. 
Again, the uncertainties are dominated by the 0.18 pixel 
uncertainties in each coordinate of the shift. We give the 
detailed fit results in Table |U 

The interval between the observations was At = 
1813.3 days = 4.96 yr, so we have a proper motion in pix- 
els (fx x ,fx y ) = (2.4+18.4, 26.3+18.4) masyr" 1 . The celes- 
tial proper motion is (/i Q , fis) = 0'.'492(— /i x , /i a ), which is 
(fJ-a^s) = (-2.4 + 18.4, 26.3 + 18.4) masyr" 1 . If we actu- 
ally want to limit the speed of SGR 1900+14, we can look 

at the distribution of the magnitude fi = ^j\i x + fiy . Our 

measurements were fi x = fi x ± a„ x arL d [i y = \± y ± y , 
although cTp.K = (Tfj,,y = dp. If we look at the magnitude 

of the proper motion [i — jl ± a^, where (i — jl x + jl 2 y , 

the probability density function for [i is: 

!M = {£) lo e-C^+AW, (3) 

where Ip(x) is the m odified Bessel function of the first 
kind (|Papoulislll99lL p. 140). From numerical integra- 
tion, we find that the mean, mode, and median are all 
close to each other. For SGR 1900+14, the mean is 
33masyr _1 (vs. fi = 26masyr _1 ), and we find a 90% 
upper limit to /i to be 54masyr~~ 1 , or a transverse veloc- 
ity v± < lSOO^kms -1 . 

4.2. IE 2259+586 

We followed the same procedure as for SGR 1900+14 
above to find the reference sources for IE 2259+586. 
Here, though, we could only identify 6 sources common 
to both epochs (a seventh source was identified, but it is 
right near the gap between detectors and the positions 
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Fig. 7. — Residual position differences between epochs 2 and 1 reference sources used to compute the transformation for SGR 1900+14 
(Left, with 27 reference sources) and IE 2259+586 (Right, with 7 reference sources); the reference sources are labeled according to Tables [2] 
and [3] The differences plotted are the residuals left after fitting for 2-parameter (i.e., shift only) transformations. The filled circles at the 
centers show the uncertainty on the resulting transformations, which in this case is the uncertainty on the mean of the residuals. 







TABLE 4 








Proper 


Motion Measurements for SGR 1900- 


-14 AND IE 2259- 


-586 


Parameter 




SGR 1900+14 




IE 2259+586 




2-paramctcr 


3-parameter 


2-parameter 


3-paramctcr 



Transform x 2 /DOF 
(Ax, Ay) (pixels) a 
(Hx,Hy) (masyr" 1 ) 13 
ft (masyr _1 ) c 
(fi) (masyr _1 ) d 
H90 (masyr" 1 ) 
(vj_) (kms- 1 ) 1 
r>±,90 (kms -1 ) 6 



37.8/52 
(0.02 + 0.19,0.26 + 0.19) 
(2 ± 18,26 ± 18) 
26 
33 
54 
780d 5 
1300d 5 



37.7/51 
(0.02 + 0.27,0.26 + 0.27) 
(2 ± 26, 26 ± 26) 
26 
39 
69 
930d 5 
1600d 5 



14.1/10 
(-0.14 ± 0.25, -0.45 ± 0.25) 
(-11 + 20, -35 + 20) 
36 
42 
65 
600d 3 
930d 3 



13.1/9 

-0.15 + 0.33, -0.39 + 0.34) 
(-11 + 26, -31 + 26) 
32 
44 
73 
630cZ 3 
1000d 3 



Shift in pixels between epochs. 
Proper motion in the (a; s ky,?/sky) frame between epochs. 



c Magnitude of the proper motion, computed from ^1 fj,"^ + (j,^ . 

^ Mean of the magnitude of the proper motion, from the distribution in Eqn. [3] 

c 90% upper limit to the magnitude of the proper motion, from the distribution in Eqn. [3] 

Transverse velocity corresponding to (p), the mean of the magnitude of the proper motion, assuming a distance of 5ds kpc (SGR 1900+14) or 3^3 kpc 
(IE 2259+586). 

s Transverse velocity corresponding to the /190, the 90% upper limit on the proper motion, assuming a distance of kpc (SGR 1900+14) or kpc 
(IE 2259+586). 



disagreed by « 3a between the two epochs). While for- 
mally enough to determine a robust transformation, this 
is fewer than we might like to have. One problem here 
is that we are using ACIS-S, which has a smaller field of 
view. Here IE 2259+586 is close to the center of a CCD 
(ACIS-S3), and we used sources from both ACIS-S3 and 
ACIS-S2. We chose to do the second epoch observation 
with ACIS-S as well, to best reproduce the first. There 
were some slight problems, though. As the roll angle 



changed by 121°, the position of ACIS-S2 changed rela- 
tive to the aim-point. Therefore, some of the reference 
sources detected in the first epoch were missed entirely 
by the detector in the second epoch, and vice versa. This 
also explains why some of the reference sources have such 
drastically different count-rates between the two epochs: 
some of them move from the front-illuminated (and hence 
less sensitive) ACIS-S2 to the back-illuminated ACIS- 
S3, although in other cases the difference in count-rate 
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1 kyr ago 



Fig. 8. — Left: Radio image of the field around SGR 1900+14 (Frail, Kulkarni, & Bloom 1999; marked with a diamond), showing 
SNR G42. 8+0.6 (to the south of SGR 1900+14). The image is a 332 MHz image from lKaplan et al.l t200% ). The arrow shows the distance 
that SGR 1900+14 has moved, starting 10 kyr ago and projecting the same time into the future using our best-fit proper motion from 
Table [4] We also show 1-, 2-, and 3-<r error circles (68%, 95%, and 99.7% confidence) on the position 10 ky r ago. Right: the same but 
plotted o n a near-infrared ( 2.2 /im) image from the Two Micron All Sky Survey (2MASS; Skrutskic et al. 2006). The position of the cluster 
found by Vrba et al. (2000) is marked with a cross. The scale of this image is reduced by a factor of ~ 10 from the radio image (the box 
on the radio image is the field shown in the near-IR image), as the proper motion vector now only shows the distance that SGR 1900+14 
has moved/will move over 1 kyr. There are scale bars in the lower left corners, and the images have north up and east to the left. 

is intrinsic, and may reflect flare stars or other variable 
sources. 

We computed 2- and 3-parameter transformations. As 
with SGR 1900+14, the 3-parameter transformation did 
not provide any improvement. With the 2-paramcter 
transformation, we find the position residuals in Figure[71 
We also repeated the jackknife and bootstrap tests (§HJ). 
The residuals in Figure [7] show an outlier, but it is rela- 
tively minor and the jackknifed transformations are con- 
sistent with the mean. 

With the transformations from the six reference 
sources, we again find only an upper limit to a displace- 
ment between the two epochs: /i < 65masyr _1 at 90% 
confidence (translating to v± < 930^ kms -1 ), with ad- 
ditional details given in Table gj Unlike SGR 1900+14, 
where the reduced \ 2 f° r the 2-parameter transformation 
was less than 1.0, here it is 1.4. This appears slightly 
worrisome, and could indicate that we have underesti- 
mated our uncertainties, which would make the inferred 
proper motion less significant. However, the majority 
of the problem is caused by a single source, #6 in Ta- 
ble [31 This source deviates by w 1.5 a, and without it 
the reduced \ 2 becomes 1.1. Including source #6, it is 
still not impossible to get such a large reduced \ 2 even 
with a proper description of the data: it should happen 
about 17% of the time, so our result is not really surpris- 
ing (also recall that our tests with the HDF-N observa- 
tions gave a reduced \ 2 of 1.33 for 1450 DOF with the 
ChaMP model). Overall, the residuals for the remaining 
sources are reasonably well behaved, and given the tests 
for SGR 1900+14 and those in §[3T] we believe that our 
uncertainties are on the whole reliable. 




1 kyr from now 



4.3. Prospects for Improvement 

Our result here is formally a non-detection of the 
proper motions of SGR 1900+14 and IE 2259+586. 



Fig. 9.— Radio image of the SNR CTB 109 at 1.4 GHz, with the 
position of IE 2259+586 marke d. The image is fro m the Canadian 
Galactic Plane Survey (CGPS; ITavlor et aUl2003l). The diamon d 
shows the current position of IE 2259+586 (Hullcman et al. 2001). 
The arrow shows the distance that IE 2259+586 has moved, start- 
ing 10 kyr ago and projecting the same time into the future using 
our best-fit proper motion from Table [4] We also show 1-, 2-, and 
3-cr error circles (68%, 95%, and 99.7% confidence) on the position 
10 kyr ago. The scale bar in the lower left indicates 6', and the 
image has north up and east to the left. 

However, the direction and magnitude of SGR 1900+14 
are suggestive, pointing as they do to a nearby SNR (sec 
§ [5]) . It would therefore be interesting to confirm or re- 
fute our measurements with newer data. Beyond this, 
one must always view 2-point proper motions with sus- 
picion, as a line can always be fit to two points, and 
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only with the inclusion of a third observation does one 
begin to gain confidence in the result. The biggest limit- 
ing factor is astrometry of th e reference sources (this is 
similar to what was found bv lMotch et al . 2007a. fj): in- 
creasing the number or accuracy of those measurements 
would be ideal, but that would have to be done for all of 
the observations. Just improving the newest one would 
have little impact. If we assume similar observations 
to those presented here, we can achieve frame match- 
ing to a precision of « 0.2 pixels between a pair of ob- 
servations. To measure a proper motion of the order of 
~ 26masyr _1 , we would need to wait > 10 yrs from 
the first observation for a statistically significant (> 3 a) 
measurement. This is only considering the first and last 
data-sets, though: by analyzing all three observations 
together we can establish a more robust reference frame 
and confirm the measurement of displacement. Regard- 
less, we would need to wait until 2010 or so before this 
will be possible. Observations at other wavelengths, pos- 
sible for IE 2259 +586 at least, should b e able to achieve 
higher precision (| Cameron et al.l 120081 . P. B. Cameron 
2008, pers. comm.). 

Another approach would be to use th e High Reso- 
lution Camera (HRC: [Murray et al.l Il997h on Chandra. 
The pixel scale is almost factor of four smaller than 
that of ACIS, so the HRC has the potential for better 
spatial resolution. However, the effective resolution is 
« 0'.'5, not that different from ACIS. It also does not 
suffer from pileup, so bright sources can be observed 
without penalty. The main drawback is that because 
of lower over all efficiency few, if any, reference sources 
are de tected: iHui fc Becker] (|2006ft and I Winkler fc Petrel 
(|2007ft found only two and three sources, respectively. 
The accuracy of our present work is limited by the regis- 
tration of the separate observations which depends on the 
reference sources. This was true even for IE 2259+586, 
where pileup lim it ed th e centroiding accuracy. Even so, 
iWinkler fc Petri Ip007ft found a significant proper mo- 
tion of 165 ±25 mas yr" 1 for RX J0822-4300 usin g three 
HRC observations (also see IHui fc Becker! 120061 for an 
earlier but consistent result). 

5. DISCUSSION AND CONCLUSIONS 

Our proper motion limit can be used to limit the age 
of SCR 1900+14/SNR G42.8+0.6, assumi ng an associa- 
tion. From the 327-MHz radio data from Kaplan et aTl 
(12002ft . we represent SNR G42.8+0.6 as a ring centered 
at (J2000): a = 19 h 06 m 59!4, 5 = +09°03'05", and with 
a radius of « 9.5' (see Fig. |HJ}. SCR 1900+14 would 
have needed > 16 kyr to traverse the separation of 15', 
using the 90% confidence upper limit on the proper mo- 
tion. This is not a remarkably large age for a ma gnet ar 
(|Gaensler et al.ll200l IWoods fc Thompsonl 12006ft . Sug- 
gestively, the marginal proper motion that we find is di- 
rected to the north, and SCR 1900+14 sits beyond the 
northern rim of SNR G42.8+0.6. However, our proper 
motion is so uncertain in both direction and magnitude 
that an asso ciation between S GR 1900+14 and the clus- 
ter found by lVrba et al.l (|2000ft is also consistent with the 
data. 

For IE 2259+586, we find an upper limit to the 
proper motion. The SNR CTB 109 (see Fig. [9]), with 
which IE 2259+586 is generally taken to be associated, 
has a somewhat complicated morphology (Koth es~et al.l 



12002ft with what may be multiple shells and seems to 
be interacting with a g iant molecular cloud to the west 
(|Tatematsu et al.l fl990) . The site of the explosion is not 
particularly well defined, given the different shells and 
the impeded expansion to the west. However, there does 
not appear to be any suggestion that IE 2259+586 is 
displaced to the south relative to the explosion center, 
as our proper motion limit might suggest. The magni- 
tude of the proper motion limit is not huge: for an age 
of 1000^3 yr, IE 2259+586 would have moved 1^3 arcmin, 
so it could have still had an origin within the current 
« 30'-diameter shell for ages up to t% ~ 10. Therefore 
our data do not greatly constrain the association between 
IE 2259+586 and CTB 109 (generally assumed to be se- 
cure) . 

Taken together, we can try to use the velocities we 
measure/limit to understand if magnetars have a ve- 
locity distribution comparable to that of the radio pul- 
sar population. The proper motion of XTE J1810— 197 
(Helf and et al.|l2b07ft and the velocity limits inferred by 
iGaensler et al.l (|2001f ) are less than 500 km s -1 . For the 
magnetars considered here, though our results are less 
constraining. This is made worse by the very uncertain 
distances to these objects. If we accept the association 
between SGR 19 00+14 a nd the massive star cluster of 
IVrba et all (|2000L also see lWachter et al J 12008ft . then we 
must also accept a distance to SGR 1900+14 of ~ 15 kpc. 
This means that our 90% upper limit on the trans- 
verse velocity is 3800 kms -1 , which encompasses the ex- 
treme range of velocities posited for the magnetars (e.g. , 
Duncan fc Thompsonlll992l iThompson fc Duncan! Il993l 
1995), but the fact of the association itself would im- 
ply a small proper motion on the order of ~ 5masyr _1 . 
For IE 2259+586, our transverse velocity limit is mod- 
er ately constraining (although w ith the larger distance 
of iDurant fc van Kerkw iik 2006 it increases to almost 
2500 kms" 1 ). 

Future astrometric observations will help answer some 
of these questions. We now understand in detail what 
the limiting factors are, and a third epoch of data around 
2010 should allow definitive measurement of the proper 
motions presented here. However, the velocity situation 
will remain more uncertain, owing to the factor of 3 un- 
certainty in the distance to SGR 1900+14 as well as un- 
known projection effects. 
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